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Abstract 

We analyze the dynamics of a quantum mechanical system in interaction with a 
reservoir when the initial state is not factorized. In the weak-coupling (van Hove) 
limit, the dynamics can be properly described in terms of a master equation, but a 
consistent application of Nakajima-Zwanzig's projection method requires that the 
reference (not necessarily equilibrium) state of the reservoir be endowed with the 
mixing property. 
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1 Introduction 



The reduced dynamics of a quantum system in contact with a reservoir is 
generally described in terms of a master equation. According to a widely ac- 
cepted lore, the physical and mathematical assumptions that are required in 
order to derive such an equation are of two types: the reservoir is much larger 
than the system (in a sense that can be made mathematically precise) and 
the coupling between them is very weak. In these limits, the system has a 
negligible influence on the reservoir and the global properties of the latter 
remain unaffected during the evolution. In turns, this enables one to assume 
that the reservoir is in an equilibrium state (which can be properly defined by 
virtue of its macroscopic features — the large number of degrees of freedom). 
One should notice that a physically clear-cut distinction between a "large" 
reservoir and one of its subsystems can only be made because of their (small) 
mutual coupling. In this respect, the two above-mentioned hypotheses are not 
only consistent, but also logically intertwined. Excellent introductions to this 
subject can be found in Refs. [1,2,3,4]. 

There is, however, another important requirement that is necessary in the 
derivation of the master equation, but that is often taken for granted: a fac- 
torized form of initial conditions is used to define the dynamics. This is a 
hypothesis of initial statistical independence that is certainly less easily justi- 
fied. The objective of this article is to investigate the evolution of the system 
with correlated initial conditions. We shall see that in such Nakajima- 
Zwanzig's projection method [2,4,5,6] requires an interesting refinement and 
a characterization of the reference state of the reservoir. 

Several authors have addressed the question of the modifications that arise 
when it is not permissible to assume initially independent system-environment 
[7,8,9,10,11,12,13,14,15,16,17,18,19,20]. In such a case, the features of the re- 
duced dynamics have interesting spin-offs in relation to other issues, such as 
the (complete) positivity of the evolution [11,12,13,14,15,16,17,18]. 

We will start off by introducing notation and setting up a general framework 
in Sec. 2. Nakajima-Zwanzig's projection method is discussed in Sec. 3 and 
our main result on the evolution when the initial state is not factorized is 
given in Sec. 4. Some of the hypotheses that are necessary in order to prove 
our results are thoroughly discussed in Sec. 5. The main theorem is proved in 
Sec. 6. In Sec. 7, we analyze the consistency of the method and we conclude 
with some remarks in Sec. 8. In this article, in order to stress the essential 
features in a transparent way, we will restrict ourselves to the main arguments 
without emphasis on mathematical rigor. A more rigorous investigation will 
be presented elsewhere. 
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Since some of the issues discussed in this article might not be familiar to ev- 
ery physicist, in particular because of the mathematical background required 
to study nonequilibrium statistical mechanics, we added some detailed expla- 
nations, translating some more advanced mathematical concepts into rather 
simple physical examples. Several details of the derivations, as well as more 
tutorial issues, are discussed in Appendices A-D. The tutorial discussions are 
often adapted to the specific problems discussed in Sees. 4 and 5. 

This is the first of a sequel of two papers. Here, the general analysis is pre- 
sented. In the following article [21], hereafter referred to as Article II, the main 
theorem on the factorization of the state, as well as the hypotheses that are 
necessary for its derivation, will be scrutinized in terms of two exactly solvable 
models, in which an oscillator is coupled to a bosonic reservoir. 



2 Framework and Notation 



We assume that the total system consists of a "large" reservoir B and a "small" 
(sub)system S, so that the total Hilbert space 7itot can be expressed as the 
tensor product of the Hilbert spaces of the reservoir TCb and of the system 

H tot =H s ®H B . (2.1) 

The Liouvillian of the total system is written as 

C = C + A£ S b, (2.2) 

where A is the coupling constant and 

Co = C s ® 1b + Is ® jCb = C s + C B (2.3) 

is the free Liouvillian, describing the free uncoupled evolutions of the system 
(Cs) and of the reservoir (Cb)- In the second equality, we omitted the tensor 
product with the unit operator and, with an abuse of notation, identified Cs 
and Cb with their dilations £s®1b and 1s®Cb, respectively. In the following, 
we will always adopt such a convention, whenever no confusion can arise. It 
follows that 

[C B ,Cs] = 0. (2.4) 

Let us write the resolution of the system Liouvillian £ s m terms of its eigen- 
projections Q m , 

Cs i ^ ' ^mQm-i ^ ' Qrn 1; QmQn ^mnQm- (^-5) 

m m 
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This resolution will be used in the following. By making explicit use of the 
Hamiltonians 

H = H + \H SB = H S + H B + \H SB (2.6) 

and noticing that, by the definition of Liouvillian, Cp = —i[H,p], the super- 
operators Q m in (2.5) can be expressed in terms of the eigenprojections Qi of 
the Hamiltonian H s as 

QmP= E QiPQi, H s = J2E i Q i . (2.7) 

Ei—Ej=ui m 

We are assuming that system S is finite, dim7is < oo, and thus has a pure 
point spectrum. 

Let p(t) be the density matrix of the total system at time t. As we anticipated, 
the initial state of the total system, p , is usually taken to be the tensor product 
of a system initial state ps and a reservoir state p B , 

Po = Ps®Pb- (2.8) 

This is an uncorrelated initial state. The reservoir state is assumed to be 
stationary (with respect to the reservoir free evolution C B ) 

C bPb = (2.9) 

and thus belongs to the eigenvalue of C B . In most applications, p B = 
Z (3 " 1 e" /3 ^ B is a thermal state at the inverse temperature (3 = (/c B T) _1 with the 
normalization constant Zp. However, we have in mind more general instances, 
such as nonequilibrium steady states [22]. The main target of the present arti- 
cle is the reconsideration of the assumption of the factorized initial condition 
(2.8). 

The system state p s (t) at time t is given by 

ps(t) =tr B p(t), p(t)=e ct p , (2.10) 

where tr B : ^(Titot) — > 7^(7is) is the partial trace over the reservoir degrees 
of freedom, a linear operator from 7i(7Y tot ), the Banach space of trace-class 
operators on the total Hilbert space 7itot 5 onto 7i(7is)- In general, unlike 
pit) = e~ lHt poe lHt , psit) is not unitarily equivalent to ps and the system 
undergoes dissipation and/or decoherence. 
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3 The Projection Method 



Consider the initial-value problem 

j t p(t)=Cp(t), p(0)=p o , (3.1) 

where the initial density operator po is not assumed here to be factorized like 
(2.8). We are interested in the reduced dynamics of the system, p s (t) given by 
(2.10), with a correlated initial state p . 

The starting point of Nakajima-Zwanzig's procedure is the introduction of the 
projection operators [2,4,5,6] 

Vp = tr B {p} <g> n B = a® Q B , Q=l-V, (3.2) 

where p G Ti(Ti, to t), & £ 1~l(Hs), and Qb £ 1i(7~(-b) is a given reference state 
of the reservoir. Note that, from the normalization condition tr B Q B = 1, it 
follows that V 2 = V and Q 2 = Q. Therefore, V is the projection onto the 
space of operators of the form a ® Qb, a subspace of Ti(Ti, tot ) isometrically 
isomorphic to Ti(Hs)- In particular, 

Vp(t) = ps(t)®n B , Qp(t) = p(t)-ps(t)®n B , (3.3) 

where we used the definition (2.10). 

Since we are interested in a correlated initial state po, which is not factorized 
like (2.8), a question arises as to which state should be taken as the reference 
state f2 B [6,9,10] and whether the naive choice f2e = tr s p is in principle 
appropriate. We shall see that the situation is much more subtle than one 
might naively expect: the reference state Qb and the reservoir dynamics cannot 
be independent, but must satisfy some important properties in order to yield 
a consistent description of the physical dynamics. One of the main subjects of 
this article will be the specification of the correct state f2 B - Furthermore, our 
analysis will show that attention should be paid to the reference state Qb even 
for a factorized initial state. This corroborates and sheds additional light on 
the rigorous conditions for the derivation of the master equation [23]. At this 
moment, we take for granted only the stationarity (2.9) of f2e with respect to 
the reservoir free dynamics, namely £b^b = 0. 

Let us project the Liouville equation (3.1) onto the two subspaces defined by 
V and Q. Notice first that 

[V, C s ] = [Q, Cs] = 0, VC B = C B V = 0, [Q, C B ] = 0. (3.4) 

The first equation is a consequence of the fact that Cs and V essentially oper- 
ate in different spaces, while the second derives from (2.9) and from the char- 
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acteristic structure of the Liouvillians, tr{£p(t)} = (a direct consequence of 
probability conservation). In addition, we require that 

VCsbV = 0, (3.5) 

which is always satisfied as long as iJsB has vanishing diagonal elements with 
respect to the reservoir degrees of freedom. 

By making use of (3.4) and (3.5), the total Liouvillian can be decomposed as 

C = VC S V + QC Q + XQC SB Q + XVC SB Q + XQC SB V. (3.6) 

Therefore, the free evolutions Cs and C B leave invariant the two subspaces 
RanP and RanQ, and all transitions are driven by the interaction £sb- 



4 The Main Theorem 



The main result of this article is the following theorem, that will be proved 
in Sec. 6: for a correlated initial state p , van Hove's "A 2 t" limit [1,24] of the 
"P-projected density operator in the system-interaction picture, 

Pl (r) = limpl A) (r) = lim e^^A>p( r /A 2 ), (4.1) 

A — >() A — >u 

is the solution of 

T 

Pi(r) =V Po + J dr'Kp^r') (4.2) 
o 

with 



/C = W dtVQ m £ SB e Cot C SB e- Cot Q m V 

m J Q 

oo 

= W dtVQ m CsBe [Ca+lWm)t CMnV 

m J 

= ~ E 'PQmCsB——^ ^-CsBQrnV, (4.3) 

or equivalently, 

^Pi{t) = /Cpi(r), Pl (0) = Vp = tr B {p } <8> ^b- (4.4) 
ar 

That is, even if the initial state po is n °t in a factorized form, but rather 
there is entanglement, or simply a classical correlation, between system S and 
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reservoir B, all correlations disappear in van Hove's limit and system S behaves 
as if the total system started from the factorized initial state in (4.4) with a 
reservoir state specified below. 



Moreover, one shows that 



lim Qp(r/\ 2 ) = lim p(r/A 2 ) - tr B {p(r/A 2 )} <g> Q B = 0, 

A — ^0 A — >0 L 



(4.5) 



which makes the dynamics consistent, for no spurious term develops in the 
master equation and no correlations can appear at later times: not only the 
initial state, but also the state at any moment t is factorized in van Hove's 
limit. This supports the validity of the assumption of the factorized state, that 
is frequently applied in literature in order to derive a master equation [2,3,4]. 
The state of system S evolves according to the master equation (4.4), while 
reservoir B remains in the state f^. 

These statements are proved under the following assumptions: 

(i) is a simple eigenvalue of the reservoir Liouvillian Cb corresponding 
to the eigenvector Qb and the remaining part of the spectrum of Cb 
is absolutely continuous [23,25] (strictly speaking, for infinitely extended 
reservoir, the spectrum of Cb can be properly defined only once the sector 
has been specified: in our case, the relevant sector is that containing the 
state f^); 

(ii) the initial (correlated) state of the total system is given in the form 



where A is a bounded superoperator (i.e., Lj's are bounded operators) 
satisfying the normalization condition tr p — 1. In other words, the initial 
state po is a bounded perturbation of the state Is <8> (and belongs to 
the sector specified by it). 

Several comments are in order. First, observe that the Liouvillian of an in- 
finitely extended system can bear a point spectrum, as in hypothesis (i) (see 
for instance Proposition 4.3.36 of Ref. [26]). 

Second, the spectral properties required in hypothesis (i) imply that Qb is 
mixing with respect to the reservoir dynamics e Cst [23,25,26,27,28], i.e., 



for any bounded (super)operators X and Y of the reservoir, where X(t) = 
e~ CBt Xe Cst and (X) p = trB{ATp}. A typical and familiar example is the ther- 
mal equilibrium state of a free boson system at a finite temperature, as ex- 
plicitly recalled in Appendix A. Among other interesting mixing states, there 




(4.6) 



(X(t)Y} nB = (Xe^Y)^ 



{X)n B (Y)n B as t -> oo 



(4.7) 
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(b) 



0', // 

> 



Fig. 1. Two typical examples in which the state of the reservoir is mixing: (a) system 
S in contact with a reservoir B in thermal equilibrium at the inverse temperature (5 
and the chemical potential fi; (b) system S in interaction with a reservoir B made 
up of two subreservoirs at different temperatures f3 and f3' and chemical potentials 
fi and //, with a steady current flowing between them (NESS). In both cases, A 
is the coupling constant: a proper application of van Hove's "A 2 i" rescaling in the 
derivation of the master equation for S requires that the reference state of reservoir 
B be mixing. 

are nonequilibrium steady states (NESS), where system B consists of several 
reservoirs at different temperatures and a steady current flows among them 
[22]. These two cases are pictorially shown in Fig. 1. Both situations are within 
the scope of our analysis. 

Third, hypothesis (i) can be relaxed. In fact, the proof will make use only 
of the mixing property (4.7). Therefore, the remaining part of the spectrum 
can contain a singular continuous part: one must simply make sure that the 
continuous spectrum is transient (a property automatically verified by the 
absolutely continuous part, due to Riemann-Lebesgue's lemma). 

Fourth, the correlation between system S and reservoir B is more transparent 
if the initial state po i n (4-6) is written as 

Po = Ps ® Pb + 5p , (4.8) 

where 

Ps = tr B p , p B =tr s po- (4.9) 

The last term Spo represents the correlation between system S and reservoir B. 
The important point is that each term is a bounded perturbation of the state 
Is <%> and thus belongs to the single sector specified by it (see Appendix B 
for a clarification and a tutorial discussion of this issue). Indeed, the factorized 
part of the initial density matrix can be written as 

p s <g> p B = A'(l s <g> tt B ) = L imn{l$ ® ^B)L\ mn , (4.10a) 

i,m,n 

Limn = y/PS <8> s((f>m\Li\(pn)s, (4.10b) 
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and 



5 Po = 6A(l s ®n B ), <SA = A-A', (4.10c) 

where {\4> n )s} is an y complete orthonormal basis of system S. Clearly, 5p is 
self-adjoint, but it is not necessarily positive-definite. 

Fifth, the states of the type (4.6) belong to a very general class: basically, if one 
deals with separable Hilbert spaces, (4.6) is of the Kraus form and therefore 
represents a completely positive map (that connects any two density matrices). 
For quantum mechanical systems (with discrete spectra), (4.6) covers all the 
possible states. In more general cases (with infinitely extended systems), since 
only bounded observables can be measured, one usually deals with a sector, 
i.e. a set of states which are "normal" with respect to some reference state (cf. 
the arguments of Sees. III. 2 and III. 3 of Ref. [27]). In such a case, any normal 
states with respect to 1s®^b can be written as (4.6) with arbitrary precision. 
(See also the arguments in Ref. [25].) 

Finally, the projection V in (3.3) must be defined in terms of Qb, which is 
mixing and "contained" in the initial state p in the sense of (4.6). This is the 
criterion for a consistent choice of the reference state f2 B . We will see in the 
next section that Nakajima-Zwanzig's projection V with the correct reference 
state Qb is nothing but the eigenprojection of the Liouvillian Cb belonging to 
the simple eigenvalue 0, as suggested by (3.4). 

Note further that a state p that refers [in the sense of (4.6)] to a coherent 
superposition of two (or more) different mixing states of an infinite reservoir 
cannot be a physical state, since it is the superposition of states belonging 
to different inequivalent sectors. (Imagine, for example, the superposition of 
states with different temperatures.) Hypothesis (ii) is therefore reasonable 
from this point of view and the states p °f the form (4.6) cover diverse phys- 
ically interesting states, ranging from canonical equilibrium states to NESS, 
as far as their f^'s are mixing. 



5 The Diagonal Projection 

Before we prove the theorem, let us observe that the projection V with the 
correct reference state f2 B is nothing but the eigenprojection belonging to the 
simple eigenvalue of the reservoir Liouvillian Cb- To this end, we first show 
how the eigenprojection of acts on a state of the total system. 

Hypothesis (i) states that for any reservoir state p B in question the spectral 
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resolution of e Cst pB reads 

e £Bi P B = n pB + J e- wt dU(u)p B , (5.1) 

where n and IT(z/) are the spectral projections of £ B belonging to its simple 
eigenvalue and to its absolutely continuous spectrum {v>}, respectively. In 
particular, 

£ B n = n £ B = o, = n . (5.2) 

By Riemann-Lebesgue's lemma, one gets from (5.1) 

e CBt PB ^ n oPB = (5.3) 

in a weak sense. The last equality follows from the condition that f2 B is the 
eigenstate belonging to the simple eigenvalue 0. 

As anticipated in the previous section, the requirement of an absolutely con- 
tinuous spectrum is not really necessary to prove (5.3). One needs only the 
mixing property (4.7). Indeed, let us consider a state of the total system be- 
longing to the same sector of Is <8> fi B , i- e - a state p of the kind (4.6). For 
any operator D — J2i A ® X i: where A^s are operators of system S and Xj's 
bounded operators of reservoir B, the mixing property (4.7) of f2 B yields 

tr{L>e £s V} = E tr B[^e £B *tr s {AA(ls®fi B )}] 

i 

^ J2 tr B {A^ B } tr{AA(l s <8> ^b)} 

i 

= tr[D(tr B {p}®n B )}, (5.4) 
and in this sense, we have 

e CBt p^tr B {p}®n B: (5.5) 

for any state p of the form (4.6). By comparing (5.3) and (5.5), it is clear that 
the eigenprojection n of the Liouvillian £ B belonging to the eigenvalue acts 
on the total system as 

(l s ®n )p = tr B {p}®fi B =Vp. (5.6) 

We have thus proved that (5.6) holds even when the spectrum {z/} in (5.1) 
contains a singular continuous part, provided that the latter be transient, i.e., 
Riemann-Lebesgue's lemma hold also for it. 

Furthermore, it is interesting to observe that the eigenprojection Is ® n is 
nothing but the "diagonal projection," that extracts the diagonal part (with 
respect to the reservoir degrees of freedom) of a density operator. To discuss 
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the diagonal components of a reservoir with a continuous spectrum, consider 
a large parameter V, corresponding to the volume, so that the Hamiltonian 
H B V ^ admits a discrete spectrum {E^}, i.e., 

HP = Y,EMPM £Pr = !> P^=P^ (5.7) 

fj, fj, 

with E^ for p,^v. Then, it is easy to define the diagonal part (with 

respect to the reservoir degrees of freedom) of a density matrix of the total 
system, p^ v \ 

Since E^'s are discrete, one easily sees that 

Km f Jdte-^ * = V, (5-9) 
o 

and thus, 



„(V)_ p (v) (V) p (v)r 

= hm 1 / Are^VVpM e ^ 

T^oo T J ^ 
M>M 

= hm 1 T Ute-< h p^e< h 

T^oo T 7 ^ 


= Km 1 Jdte&'pM. (5.10) 
o 

We now define the diagonal part of p via 

1 r 

p D = lim p D v) = lim — / dte CBt p, (5.11) 



assuming that the two limits can be interchanged. Then, under the assumption 
of the ergodicity of Qb [which follows from the mixing of Qb' see Eq. (A. 4)], 
one has, for any state p of the type (4.6), 

Pd = tr B {p} <g> n B , (5.12) 

which is the right-hand side of (5.6). Indeed, for any bounded operator D = 
Y,iAi ® Xi considered in (5.4), 
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= tr B {X^ B } tr{AA(ls ® ^b)} 



tr[D(tr B {p}®n B )}, 



(5.13) 



which proves (5.12), and the eigenprojection of Cb belonging to the discrete 
eigenvalue is the diagonal projection. Notice that the ergodicity of fi B is suf- 
ficient to show that Is <8>n is the diagonal projection. A direct demonstration 
of the diagonal projection (5.12) in a simple model is given in Appendix C. 

It is now clear that the eigenprojection ls®n in (5.6) acts like the projection 
V defined in (3.3), provided the state p(t) is of the form A t (l s <S> ^b), with 
a bounded superoperator A t . Since p(t) has evolved from the initial state 
(4.6), it always satisfies this criterion. As we will discuss in the following, the 
eigenprojection l s <g> Il enables us to deal with the point spectrum and yields 
the right choice for the projection V in order to derive a master equation in 
van Hove's limit. 

Summarizing, the initial state p "contains," in the sense of (4.6) and (5.5), 
the mixing state Qb- The theorem stated in Sec. 4 will be proved in Sec. 6 
with the projection operator 



which is the eigenprojection of Is <8> £b belonging to the unique simple eigen- 
value 0. In the following, with an abuse of notation, l s <8> il will be written 
simply as II , as far as no confusion can arise. 

Some additional comments are in order. If we assume only the ergodicity of 
CIb, rather than mixing, there can be other eigenvalues different from 0. The 
point spectrum other than is out of control and, as we shall see in the 
following section, provokes the appearance of a secular term. 

It is worth noting that the mixing property of f2 B is crucial even for a factorized 
initial state like (2.8), although this point is usually not thoroughly discussed. 
(This is an interesting byproduct of our analysis, that is rather motivated by 
the study of nonfactorized initial states.) As we shall see, a wrong projection 
V would give rise to a divergence that has nothing to do with the initial 
correlation. 

In the standard derivations of the master equation for a factorized initial state, 
the projection (3.3) is defined in terms of the same Qb which is contained in the 
factorized initial state (2.8) and is usually the canonical state at temperature 
T, that is a mixing state. The criteria listed in Sec. 4 are thus satisfied and the 



^ = i s ®n , 



(5.14) 
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projection is properly chosen to be V = lg£g>ITo, provoking no problem. In this 
standard situation, the choice of the reference state is obvious: probably this 
often induces one to assume that Nakajima-Zwanzig's reference state is simply 
the canonical one (and therefore need not be characterized). However, in more 
articulated situations, such as those of NESS (see Fig. 1), the reference state 
must be chosen with care: our analysis shows that a correct reference state 
must in general be mixing. 



6 Proof of the Theorem 

Let us now prove the theorem stated in Sec. 4. By projecting the Liouville 
equation (3.1) onto the two subspaces defined by V and Q, and using (3.6), 
one gets 



-Vp = C s Vp + XVCsb Qp, (6.1a) 

j f Qp = jC' Qp + \QC SB Vp, (6.1b) 
respectively, where 

C' = C + XQC SB Q. (6.2) 



By formally integrating out the second equation and plugging the result into 
the first one, one gets the following exact equation for the P-projected operator 
in the interaction picture [6] 

, t 
j t e- Cst Vp(t) = \ 2 J dt'e- c ^VC S Be c ^C SB Vp(t') + Xe'^VC^Qp,. 
o 

(6.3) 

The last term in the right-hand side represents the contribution arising from a 
possible initial correlation between system S and reservoir B. We will show that 
this term dies out in the weak-coupling limit A — > with the scaled time r = 
A 2 t(> 0) fixed, provided the projection V is chosen to be the eigenprojection 
Is <S> n belonging to the simple eigenvalue of C B . To this end, consider the 
density operator 

p[ X \r) = e' c ^ x2 Vp(r/\ 2 ), (6.4) 
introduced in (4.1), that for any nonvanishing A satisfies 



^-p ( l X \r)= [ dte- c ^ V c SB e^/^)c SB Ve^p?\\h) 
dr J 
o 

+ \e- c ^ x2 VC SB e c 'o^ 2 Q P(h (6.5) 
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with the initial condition 

p{ X \0)=V P0 . (6.6) 

By integrating (6.5), one gets 



r'/A 2 



p W {T)=Vpo + J dT ' J dte-^'/ x2 VC SB e c ^'^^C SB Ve^p{ X) (X 2 t) 
o o 

+ I J dr> e-^'/ x2 VC SB e^'/' 2 Qp . (6.7) 
o 

The second term is rearranged as 



l l\2 



Jdr' J dte-^'^VCsBe^'^CsBVe^p^iXH) 



r/A 2 \H 

= J dtj dr' e- c ^VC SB e c ^'/^C SB Ve c ^/ x2 p{ X \r') 
o o 

r/A 2 

= Jdr' J dte-^VCsBe^'^CsBVe^'^p^ir') 

r'/A 2 

(r-r')/A 2 

= Jdr' J dt e-^^'/^VCsBe^CsBVe^'/^p^ (r') 



T 

= E / rfr'/— ^'/^W(r - r')p[ X \r'), (6.8) 

m,n q 

with the memory kernel 

r/A 2 

K^(t)= / dtVQ m £ S BQe {c » +tuJm)t £sBQnV. (6.9) 



The last term in Eq. (6.7), which is relevant to the initial correlation, reads 

lW(r) = \ Jdr'e~ c ^/ x2 VCsBe c " T ' /x2 Qpo 
o 

r/A 2 

= \J2 [ dtVQ m £sBe {c '» +tuJm)t Qpo- (6.10) 
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In conclusion, 

T 

Pi (A) (r) = V Po + £ / dr'e^-^'^^Kr - r')p[ X) (r') + 2<*> (r). (6.11) 



We can now analyze the van Hove limits of the memory kernel /C^(r) and of 
the initial correlation X^ x \t). Both limits can be computed at the same time 
if we consider the van Hove limit of the operator 

H$(j)= J dt Qe {c '" +iu)m) \ (6.12) 
o 

since both /C^(r) and l^ x \r) contain it. The analysis of the properties of 
1Z!$(t) is given in Appendix D, where it is shown that, irrespectively of the 
point spectrum of C , one obtains 

^ = lim^g)(r) = - Q (r>0), (6.13) 

provided that V = Is <8> n . The expression (6.13) is our key formula: it imme- 
diately leads us to the conclusion that, in van Hove's limit, the memory kernel 
^mn( r ) i n (6-9) is reduced to a Markovian generator 



^ /C£ = -VQ m C SB ^— — ^CsbQuV, (6.14) 

Lq + lLO m - 0+ 

while the correlation term X^(r) in (6.10) disappears 



JM ( T ) = A£ VQ m C SB n^ (r) Qp ^ 0, (6.15) 

rrt 

so that the reduced dynamics (6.11) becomes 

T 

p[ X \t)^V Po + Y: I dr' 5 mn K^ nPl {r') 

m,n q 

T 

= V Po + j dr'Kp^r'). (6.16) 
o 

The master equation in van Hove's limit therefore reads 
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— pi(r) = /Cpi(r), 



(6.17a) 



Q 



(6.17b) 



rn 



rn 



C + iuj m - 0+ 



We have thus proved (4.2)-(4.4), the first part of the theorem in Sec. 4. In 
van Hove's limit, the density matrix evolves as if it started from the initial 
condition Vpo, even when Qp ^ 0: the initial correlation dies out immediately 
(at t = + ) in the scaled time r. The contribution originating from the initial 
correlation between the system and the reservoir disappears in the scaling 
limit and one is allowed to start from an initial density matrix in the factorized 
form (4.4). From a physical point of view, the factorization Ansatz described 
above simply means that the "initial" correlations between the system and its 
environment are "forgotten" on a time scale of order A 2 in r (which is very 
small when compared to the timescale of the evolution of the system). We 
shall see in Article II [21], by looking at a specific example, that the problem 
of the relevant timescales must be tackled with it also involves locality 

issues related to the "size" of the local observables of the reservoir. 

It is important to note that, if we choose a different projection V from the 
eigenprojection of Cb belonging to its vanishing eigenvalue, and fail to appro- 
priately remove the point spectrum of Cb, Eq. (6.13) does not hold anymore, 
and H$(t) diverges in van Hove's limit. [See, for example, Eq. (D.6) in Ap- 
pendix D.] /C^(t) and X^{t) accordingly diverge and the van Hove limit of 
the master equation does not exist. This is because a "wrong" projection V 
would project the reservoir onto a wrong (in general non-stationary) state, so 
that the system evolution would develop a secular term in r/A 2 [like in Eq. 
(D.6)]. It is remarkable that such a secular term appears as a consequence of a 
sloppy application of Nakajima-Zwanzig's projection method: in a sense, the 
very method makes sense only if applied to the "correct" mixing state. 

As stressed at the end of the previous section, the reference state of the reser- 
voir must be mixing (and not, e.g., simply ergodic), in order that no discrete 
eigenvalue different from exists. Otherwise, the point spectrum (except 0) is 
out of control and again the emergence of secular terms is inevitable. More- 
over, the mixing property of the reservoir is crucial even for a factorized initial 
state like (2.8). Indeed, careful scrutiny of the proof shown in this section shows 
that the divergence of /C^(r) has nothing to do with the initial correlation. 
In this sense, although our original motivation was the study of the factor- 
ization Ansatz, the results have more general validity. The reason why this 
was unnoticed so far is the following. In the standard textbook derivations 
of the master equation, a factorized initial state is assumed and the projec- 
tion (3.3) is defined in terms of the same reference state Qb that is contained 
in the factorized initial state (2.8): in practically all examples, this is taken 
to be the canonical state at temperature T, which is clearly mixing, and all 
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the criteria listed in Sec. 4 are satisfied. However, recently, more complicated 
physical situations are attracting increasing interest, such as those related to 
NESS shown in Fig. 1(b). These cases call for a characterization of the state 
of the reservoir, that is mixing and satisfies the hypotheses in Sec. 4 required 
to prove the theorem [22]. 



7 Factorization at All Times 

Time t — has no particular status: one can prove the same master equation 
(4.4) with a factorized initial state for any "initial" time to- This means that 
the limiting dynamics is such that the density matrix remains factorized at all 
times: the reservoir state does not evolve, while the system state follows the 
master equation (4.4). In order to show this, one must prove the validity of 
Eq. (4.5) in van Hove's limit. Then, clearly, no spurious term will develop in 
the master equation and no correlations can appear. 

By integrating (6.1a) and inserting it into (6.1b), we get the following integral 
equation for Qp(t), 

t 

Qp(t)=e c '° t Q Po + \ J dt'e c o^QjC SB e Cst 'Vp 
o 

t f 

+ X 2 J dt'J dt"e c '^QC SB e c ^'-^VC SB Qp(t"), (7.1) 
o o 

which is rearranged, by interchanging the integrations in the last term, to yield 
in the scaled time r = \H 

Qp(r/\ 2 ) = Qe c '° T / x2 ( Qp - \'£K£\-T)£ SB Q m Vp 

~ E IdT'e-^'^TZ^ir' - r)CsBQ m VC SB Qp(r'/\ 2 )) , (7.2) 
m o / 

where 1Z$(t) is defined in (6.12). Under the proper choice of the projection 
V, the kernel TZ$(t) is bounded for any r and A (even for A — > 0), as already 
discussed in the previous section and in Appendix D, and it is possible to show 
that (see Appendix D) 

lim Qe c '« T,x2 = 0. (7.3) 

Then, the integrand in the last term in (7.2) is bounded, all the terms in the 
parentheses are finite, and the prefactor in (7.2) vanishes according to (7.3), 
which proves (4.5), the second part of the theorem in Sec. 4. 
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It is worth noting that the interaction between system S and reservoir B is 
not essential to the factorization; the free evolution eliminates the correlation, 
and the reservoir relaxes into the mixing state Q-b- Indeed, for any state po of 
the total system of the type (4.6), we have 

e c t po e £st nopo = e c s t trB{po} ^ 

where the contribution of the absolutely continuous spectrum decays out due 
to Riemann-Lebesgue's lemma. [See Eqs. (5.1) and (5.5).] In the rescaled time 
r, the factorization is very rapid, and the total system looks factorized at any 
moment (if the observables that one can measure on the reservoir are local 
enough: a concrete example will be discussed in the following Article II [21]). 
Summarizing, reservoir B relaxes into the mixing state through its own 
free evolution, yielding the factorization of the state of the total system, while 
system S dissipates through the interaction: a remarkable and consistent global 
view. 

It is also interesting to compare the present result with Bogoliubov's view 
on the classical gas dynamics [2,29,30]. According to this view, "molecular 
chaos" erases a large amount of initial information and the system reaches the 
so-called kinetic stage, where the one-body distribution function governs the 
evolution of the whole system and obeys the Boltzmann equation, irrespec- 
tively of the initial conditions. In the present case, the initial loss of system- 
reservoir correlations corresponds to the information loss due to molecular 
chaos, and the stage described by the master equation in van Hove's limit 
corresponds to the kinetic stage. 



8 Concluding Remarks 

We analyzed the assumption of factorization of the initial state in the dynam- 
ics of a quantum mechanical system in interaction with a reservoir. In van 
Hove's limit, the dynamics can be consistently described in terms of a master 
equation, but a correct application of Nakajima-Zwanzig's projection method 
requires that the reference state of the reservoir, in terms of which the projec- 
tion operator is defined, be endowed with the mixing property. If the reservoir 
dynamics is not mixing, the evolution develops secular terms. In the above 
discussion, one implicitly assumes that the van Hove limit (4.1) exists and is 
the solution of Eqs. (4.2)-(4.4). 

In the standard situation, when a small system is coupled to a reservoir at a 
given temperature, the choice of the reference state is obvious and is simply 
the canonical state (that is also mixing). In more articulated situations, such 
as those of NESS, the choice of the reference state requires care and need to 
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be characterized: our analysis shows that a correct reference state must in 
general be mixing. 

This Ansatz yields the standard procedure and the usual master equation. As 
a byproduct, one observes that the mixing property is crucial even when the 
initial state is assumed to be factorized, otherwise the presence of a secular 
terms is inevitable. In this respect our analysis, that was originally motivated 
by the study of the assumption of initial state factorization, has more general 
validity. 

We shall see in the following Article II [21], by close scrutiny of some explicit 
examples, that Markovianity becomes a valid approximation for timescales 
that depend both on the form factors of the interaction and on the spatial 
extension of the local observables that can be measured on the reservoir. This 
will corroborate and complement the general findings discussed in this article. 
Other issues, such as the spin-offs on the complete positivity of the dynam- 
ics [11,12,13,14,15,16,17,18], the mathematical conditions at the origin of van 
Hove's limit [31], as well as the consequences of higher order corrections [32] 
and their interplay with initial correlations will be investigated in a future 
paper. 
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A Mixing Property 



The free Liouvillian of an infinite reservoir £ B has a point spectrum [23,26,27]. 
The clue for our problem is to handle it properly, by making use of the right 
projection V. The mixing property of the reservoir plays an important role in 
this context. Let us hence briefly recall these notions. 

A.l Mixing Property and Spectrum of the Liouvillian 

The state VLb is said to be mixing with respect to the reservoir dynamics e Cst , 
if the correlation function of any two bounded operators of the reservoir, X 
and Y, behaves as (4.7) [23,25,26,27,28]. The ordinary canonical equilibrium 
state of free bosons at a finite temperature is a typical mixing state. Other 
interesting examples are nonequilibrium steady states (NESS). 

It is important to observe that the mixing property is strongly related to 
the spectral properties of the Liouvillian £ B [28]. Let us consider a reservoir 
state pb = Ab^b (tre Pb = 1) related to a mixing state Qb by a bounded 
superoperator A B , in the sense of (4.6). Then, by setting Ff2 B = A b Qb = Pb 
in (4.7), the mixing property of Qb translates into 

(X(t)) PB -> (X) Qb as t^oo, (A.l) 

and in this sense, 

e CBt p B -> as t -> oo, (A.2) 
i.e., the state pb relaxes towards the mixing state Qb- 

Let us consider the spectral resolution of e £l3 'pB, 

e £B V B = E e' iVj %PB + / e- M dU(u)p B , (A.3) 

where IT,- and n(z/) are the spectral projections of Cb belonging to its possi- 
ble discrete eigenvalues {uj} and to its absolutely continuous spectrum {u}, 
respectively. The second term, representing the contribution of the absolutely 
continuous spectrum, decays out as t — > oo due to Riemann-Lebesgue's 
lemma, but the first term, the contribution of the point spectrum, survives. 
Property (A.2) means that (within the class of states of the form p B = A b Qb) 
there exists the only simple eigenvalue of the Liouvillian £ B with the eigen- 
projection n satisfying n pB = ^b : mixing forbids the existence of other 
eigenvalues than 0, reducing (A.3) to (5.1), and the eigenvalue is not degen- 
erated within the single sector specified by Qb- Furthermore, since U(is)Qb = 
[remember the orthogonality Hon(z/) = n(z/)Ilo = 0], the mixing state Qb is 
a stationary state with respect to the reservoir dynamics e Cst , i.e. Eq. (2.9). 
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If we only require that be a simple eigenvalue of Cb, but we do not care 
about the rest of the spectrum, then property (4.7) only holds in a Cesaro 
sense, namely 

1 T 

hm - J dt(X(t)Y)a B = (X)n B (Y)n B , (A.4) 
o 

and f2 B is called an ergodic state. Ergodicity suffices to show the stationar- 
ity (2.9) of f2 B ; the mixing property is not necessary [23,25,26,27,28]. The 
uniqueness of the point spectrum 0, however, is equivalent to weak mixing 
[23,26,27,28] / T dt | (X(t)Y)n B - (X)n B (Y)n B \/T -> and this implies that 
the remaining part of the spectrum can also contain singular continuous com- 
ponents. In such a case, one might conjecture that the theorem of Sec. 4 is 
still valid, but in a weaker sense. 



A. 2 A Solvable Example 



Let us demonstrate the mixing property (4.7) in a solvable example. Let us 
consider the reservoir Hamiltonian 

H B = J dojub% (A.5) 
and the reservoir dynamics starting from the initial state 

p B = A B pw = J ' duj J dJ w^'bipwbw* (A.6) 



with 

Pw 



-^-exp(- J du J fa' blWuu/bS) , (A.7) 

where b^ (&£) is the annihilation (creation) operator of the reservoir, satisfy- 
ing the canonical commutation relation [^fyy] = S(u — uu'). The states p B 
and p w are normalized, tr B p B = 1 and tr B pw — 1> with the normalization 
constant Z w . w^u' is a bounded and Hermitian positive matrix (w uu > = w^, u ), 
and the Gaussian state pw is perturbed by the bounded superoperator A B . 
Furthermore, W WUJ > is Hermitian (W^/ = W*/ w ) and consists of W<^,, that is 
proportional to 5{uj — a;'), and the remaining part VV<^y, 

= + HLy, = W(u)6(u - J). (A.8) 



The Gaussian state pw is fully characterized by the two-point function 

Nww = (&1/Upw = N(u)6(u - u') + AL,,, (A.9) 
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where the first term, proportional to 5{uo — uo'), is the expectation value of the 
number operator in the state pw , 

•A/2 = <^>pw = N(u)S(u - a/), JV( W ) = eW{ ^ _ - (A.10) 

with 

Pw = ^ exp (~/ ^O^M 6 ") > tr B pw = 1- (A.11) 
In fact, pvv is different from pw by a bounded operator L B , as 

p w = LbPwo^L ( a - 12 ) 

where 

1/2 



= P^Pwf = V Texp ( - / d(3 bh'^We^b ) , (A.13) 



o 



with T denoting the anti-chronologically ordered product and 

b^e-^We^b = J duo J dw'ble-^W^e^'K, (A.14) 
Hence, by applying Wick's theorem, the two-point function (A. 9) reads 



ACo/ = (L^IA^b)^ = (^^)pw (4^b) PW() + • • • • (A. 15) 

The first term is A/^ in (A. 10), since (L b L b ) PWo = tr B pw — 1, and the other 
terms, defining A/^', are bounded functions, not proportional to 5{uj — a/). 

Let us now take any two bounded operators of the reservoir, of the form 

X = Jduj daj'blX uu ,b u ,, Y = jdcoj du'bly^b^, (A.16) 

and observe how the mixing property (4.7) emerges. In this case, the correla- 
tion function reads 



(X(t)Y) PB = (e lH ^Xe- iHBt Y) 



PB 



= J dco 1 J dco 2 J doo 3 J doo A J duo J duo' A^^X. 

x (b^lb^bl^X)^ 1 -" 2 *- (A.17) 

By applying Wick's theorem, {b^'b^b^bl b^b^) Pw can be expressed in terms 
of two-point functions (A. 9), and by applying Riemann-Lebesgue's lemma, 
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only the contribution of 8{uj\ — 0U2) in (^^lew survives in the long-time 
limit, to yield 

(X(t)Y) 

p B > f duii J did 3 J did 4 J did J duj X^^y^^w^i 

Pw 

= (X) PW() (Y) PB . (A. 18) 

Therefore, if pe = Pw ( Le -> C u < = without the perturbation A B ), Eq. 
(A. 18) is nothing but the definition of mixing in (4.7), and pw given in (A. 11) 
is an example of mixing state. The canonical state pw with W(uj) = /3u, is 
thus a typical mixing state. We have demonstrated (A. 18) with the specific 
observables X and Y in (A. 16), but this example helps us understand how 
mixing emerges for general observables. 

It is important to note that, in this appendix, we have considered only the 
reservoir dynamics generated by the reservoir Hamiltonian Hb, without any 
interaction. The free evolution is responsible for the mixing. 

The mixing property (4.7) is demonstrated here with the thermal equilibrium 
state (A. 11). It is also possible to prove it for the NESS depicted in Fig. 1(b). 
See Ref. [22] for details. 



B Inequivalent Sectors 



The main purpose of this tutorial appendix is to clarify that different mixing 
states belong to different sectors which are inequivalent to each other, and any 
bounded perturbation on a mixing state does not bring it to a different sector. 
Let us demonstrate the inequivalence of the sectors with an explicit example 
that captures the essence of the inequivalent representation. 

In order to analyze an infinitely extended system, let us begin with a free 
bosonic gas in a ID box whose size is specified by a parameter £, and then 
take the continuum limit I — > 00. We consider two canonical states flp and Qp> 
with different temperatures and compute the overlap between them through 

1 /2 1 /2 

the quantity tiBx^p f^/, }. Note that the canonical state is a mixing state as 
shown in Appendix A. 

In the finite box, momentum k is discrete, and the Hamiltonian of the bosonic 
gas is given by 

H B = ^fbf, (B.l) 

k 
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where bf 1 and satisfy the canonical commutation relation [b^\ 
The canonical state at the inverse temperature (5 is given by 

Z /3 fc 

and the overlap between Vtp and fi^/ reads 



(B.2) 



tr B {fif^ 2 } = n 



^(1 - e -^*)(l-e-^*) 



fc 

= exp 



1 _ e -/3w fe 

-- Y In { - - } - 

2^ (l-e-^)(l-e-^) 



(B.3) 



where f3 = (f3 -\- f3')/2. The exponent in (B.3) is easily shown to be less than 
zero for (3 ^ f3' and equal to zero for f3 = f3' . In the continuum limit i — > oo, the 
summation in the exponent is replaced with an integral as J2k ~^ (^/2vr) / dk, 
and the overlap is reduced to 



tr B {fi,!/ 2 fi«/ 2 } = exp 




dk In 



'1 _ e -/3^fc)2 



1 _ e -/3w*)(l - e ~/3'^) 

as £ — > oo, 



(B.4) 



which means that canonical states with different temperatures do not overlap 
and belong to inequivalent sectors. 

Any bounded perturbation does not change the situation: a state tip, which is 
different from a canonical state fig only by a bounded superoperator, belongs 
to the same sector as that of the canonical state flp and does not overlap with 
a canonical state Qpr with different temperature. Consider, for example, a 
state tip, whose square root fi^ 2 is different from fi^ 2 by a bounded operator 
Kb (or Lb) as 



tif = K B ttfK B 



L-B&p 2 



fiJ /2 4, L B = K B Q /2 K B Q: 



(B.5) 



with the normalization conditions tr B rip — 1 and tr B Qp = 1. We again begin 
with a finite I, so that the overlap between fi^ and Qpr now reads 



i o i -> , > I - -r w(l-e-^)(l-e-^) 

tr B {fiJ /2 fi^ 2 } = tr B {LBfiJ /2 fi^} = II ^ r^ZS- L (L 

k 1 ~~ e 



(B.6) 
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where (Lb) 3 is the expectation value of Lb in the canonical state at tempera- 
ture (3. Note that (Lb) 3 is finite even in the continuum limit i — > 00, since Lb 
is a bounded operator. Therefore, exactly the same argument as (B.4) applies 
to this case and leads to the conclusion 

tr B {jy n^, } -> I o _^ ^ as £ -> 00, (B.7) 

i.e., Q/3 does not overlap with Qpr and belongs to the sector equivalent to Qp. 



C Diagonal Projection 



Let us confirm the property of the diagonal projection in (5.12) with an explicit 
example. We consider the same model as in Appendix B and observe how the 
diagonal projection acts on a reservoir state tip = Ab^/3, which is different 
from the canonical state Qp at the inverse temperature f3 only by a bounded 
superoperator A B , and therefore belongs to the sector equivalent to Vtp. 

We begin with a ID bosonic gas in a finite box, and then take the continuum 
limit t — > 00. For a finite £, the diagonal projection P D is defined as (5.8) and 
is given in this case by 

Pbfy, = £ IK})(K}|^|K»(K}|, (c.i) 

{rik} 

where ({n^}) = \n kl nk 2 . . .), n k being the occupation number in mode u) k . To 
be explicit, let us take a reservoir state 

n g = A B a, = £ «,» »f> V» = £ «,» ^^n, (c.2) 

with a Hermitian (wj^l = w[% ) positive matrix, and consider the expectation 
value of a bounded operator of the reservoir 



k,k' 

in the projected state V^VLp: 



y = E^«> (C3) 
k,k' 



ti B {YV D ^}=Y / ({nk}\Y\{n k })({n k }\n f3 \{n k }) 

{«fe} 

= E (K}|^IK})(K}|^b|K})(K}|q^|K}), (c.4) 

{n fc } 
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where L B = J2k,k> w kk' e ^ k '^ b k ) ls the bounded operator acting on the left 
side of ftp in (C.2). Note that Qp is diagonal with respect to the basis 
and the diagonal elements are given by 

({n k }\Qp\{n k }) = l-e-f*^™ (C.5) 

with Zp given in (B.2). The diagonal elements of Y and L B are easily evaluated 
to be 



({n k }\Y\{n k }} =52yjgn k , ({n k }\L B \{n k }} = £^e^n fc , (C.6) 

k k 

respectively, and Eq. (C.4) is 

tr B {YV D n p } = ^ E T,yki^i%^'n k n k ,e-^'' n ^''. (C.7) 

Z /3 {n k }k,k> 

Notice here that 

M = i- £ n^e-^— » = (<"^<"^ <* * l\ (C.8) 
where 

(n k )p = {bfbf) p = (nl) fi = {{bfbff) p = 2{n k )l + (n k )p, 

(C.9) 

so that Eq. (C.7) is decomposed into two terms, 

tre^} = E )S4^(4(^ + E^(4 (CIO) 

We are now in a position to take the continuum limit £ — > oo by recalling the 
correspondence 

E ^ ^ / ^ ^ \Zy 6fc ' y *2 ^ y^**'' ^ y w "'- ( an ) 

where and satisfy the canonical commutation relation [b k , b\,} = 5(k — k'), 
and y kk > and w kk i are assumed to be bounded functions. The relevant quantity 
now becomes 

ti B {YV D n p } = JdkJ dk'y kk w klkl e^'{n k ) p {n k/ )p 

2ix 



+ 



Z7T f 

-jj dky kk w kk e^(nl)p (C12) 
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for large £, and the second term disappears in the continuum limit £ — > oo. 
We finally obtain 

tr B {YVv(lp}^(Y) p (L B )p = (Y) p tr B np as £ ^ oo, (C.13) 
which yields the formula for the diagonal projection, 

PoClp = n p . (C.14) 



The action of the diagonal projection Vn on the total system is now readily 
understood. Consider, for example, a state of the total system 

p = A(l s ® fi,?) = ]T j ' dk J dk'w ijM/ Sibl(as <g> fy^Sj, (C.15) 

hi 

where as is any positive operator of system S, S^'s are system operators, and A 
is a bounded superoperator acting on Is <S>flp. Take an operator of the system, 
A, and an operator of the reservoir, Y given in (C.3). Starting with a finite £, 
the expectation value of the operator D = A <g> Y in the projected state P D p 
is 



ti{DV D p}=J2(M\Y\{n k }}({n k }\tT S {Ap}\{n k }} 

{n-k} 

= EE trsl^as^j}^^,^^' Kn^, (C.16) 

i,j k,k' 

which is reduced, in the continuum limit £ — * oo, to 



tr{DV D p} j dk j dk ' tos{AS i (JsSj}ykkWij,k'k'e l3Uk ' (n k )p(n k <)p 
hi 

= tr s [Atr B {A(l s ®^)}](r) /3 

= tr[ J D(tr B {p}®^)], (C.17) 
reproducing (5.12). 



D Key Formulas for the Theorem 



Here we prove the key formulas (6.13) (with its counterpart for r < 0, which is 
necessary in Sec. 7) and (7.3), and see how the proper choice of the projection 
V is crucial. 
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D.I The van Hove Limit of Kg* (t) 



Let us analyze the kernel K$(t), defined in (6.12). In this appendix, the 
eigenprojection l s <S> II is written simply as n . 

We start by noting that in van Hove's limit, for n > 2, 

r/A 2 

\ n Kg(r) = A™ J dt Qe^+'^t as A -> (n > 2), (D.l) 

o 

irrespectively of the spectrum of C' Q . 

Second, the following observation will be important: the convolution 

t 

J dt'e^-^QCsBQe^' (D.2) 
o 

is bounded for any t, provided the point spectrum of Cq is removed by the 
projection Q = 1 — V with V = n . Let us look at the Laplace transform of 
this convolution (for t > 0), 

1 -QCsbQ—^t- (D.3) 



s - C s - C' 

Neither l/(s — £ ) nor l/(s — £q) has a singularity on the right half plane 
Res > 0. If Co and £q have common eigenvalues (along the imaginary axis 
Res = 0) that are not projected out, these would give second order poles and 
yield linearly diverging functions of t (for large t) after the inverse Laplace 
transform; otherwise, the convolution decays or just oscillates. On the other 
hand, if the point spectrum of C is removed by the projection Q, such a 
coincidence between the point spectra does not happen and the convolution 
(D.2) is bounded for t — ■> oo, irrespectively of the point spectrum of C . 

Now by using 

t 

+ \J dt'e^-^QCsvQe^', (D.4) 



we expand the relevant quantity (6.12) as 

r/A 2 r/A 2 



Kg{T)= J dt Q e ( C 0+iu>m)t + A J dt J dt > e {C +iu m )t e -C t> Q£ SB Q e ^t' . 


(D.5) 
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The first term is decomposed into two parts by the projections Ho and Il c = 
1 — n , and the integrations are easily carried out to give (for r > 0) 



r/A 2 r/A 2 r/A 2 



J dt Qe {Co+iuJm)t = J dt QU e {Cs+iuJm)t + J dt Qn c e {Co+iuJm)t 

e (C +iui m )T / \ 2 _ J 





r/A 2 



]T f dt Qn Q n e* Um - Un '> t + QU 



C + iuj m — 0^ 

!(w m -w„)r/A 2 _ y 



(C +iu>m)T I A 2 _ 1 

£ + lUm - 0+ 



which shows that the only possible divergence of the relevant operator in 
(6.12) in van Hove's limit A — > stems from the point spectrum of Cq (i.e. the 
first term of the last expression). This divergence results in the divergences of 
both the memory kernel /C^(r) and the initial correlation X^(r). However, 
ifV — Il , the divergent term disappears due to QU = 0, and we have (for 

T > 0) 



/ \2 

/JC +lU) m )T / \ 2 _ 1 Q 
dt = Q— — 7^^-^—^ as A ^°> 

Lo + %oj m - 0+ Lq + iu m - 0+ 

(D.7) 

by noting the formula 



e (C +iU!m)t 

lim n c - = 0, (D.8) 



which is valid in the sense of distributions. The second term in (D.5) can be 
manipulated to yield 
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r/A 2 t 

A J dt J dt' e ^o+iu>m)t e -C t' Q£ SB Q e C[/ 



r/A 2 r/A 2 

= A J dt' J dte^+^e-^'QCsBQe^ 1 ' 



o v 



Q 



= A- e^ri* f dte £«(r/xi-t) QCsBQe c> t 

C + lU) m — 0+ J 



o T/x2 

- — — — -£ SB / rftQe^"^)'. (D.9) 



The integral in the first term is the convolution (D.2). Since the point spectrum 
of Co is removed by the projection Q, this convolution is bounded for r/A 2 — > 
oo. 

In summary, with the choice of the projection V = n , Eq. (D.5) is arranged 
into the recurrence formula 



e (C () +ico m )r/X 2 _ 1 

C + iuj m - 0+ 

r/A 2 



+ A- e luJmT/x2 I dte^^'^QCsBQe^ 

C + %oj m — 0+ J 

-A-^— — £sB^ } ( r ) (r>0), (D.10) 

L Q + iu m - 0+ 



where the first term converges to (D.7) and the second term vanishes in van 
Hove's limit. Therefore, by iterating the above expansion twice, we arrive at 



Q ( Q_ 

\Co + itO m 

r>0), (D.ll) 



lim7^)(r) = -— - ( Q £sB) 3 lim A3^)( r ) 

A^O Lo + lOJm — U + \Lo + lUJ m — + / A^0 

Q 



Co + iu m - 0+ 
which is Eq. (6.13) of the text. 

A similar argument applies to the case r < to yield 
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e (£ +ioj m )r/A 2 _ y 

C + iuj m + 0+ 

Q 



+ A 
- A 



C + iuj m + 0+ 
Q 

C + iuj m + 0+ 





Affi^gV) (r<0), 



and 



Q 



Co + iuj m + 0+ 



(r < 0). 



(D.12) 



(D.13) 



D.2 The van Hove Limit of Qe £ o T A 2 



Using the expansion (D.4), we have 

Qe C> T/X* = Qe C T/X* + X I dte ^(r/X 2 -t)Q CsB Q e C' t^ (D14) 



The convolution in the second term is the same as that discussed in (D.2), 
which is bounded for any r and A, provided the projection V is the eigenpro- 
jection Il , and hence the second term vanishes in van Hove's limit A — > 0. 
Since the right projection Q — 1 — n removes the point spectrum of C , 
the first term disappears as r/A 2 — > oo due to Riemann-Lebesgue's lemma. 
Therefore, Qe c » T l x decays in van Hove's limit, yielding (7.3). 
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